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ABSTRACT. There are many ways of programming an actuator or effector 
for movement between th same two points. In the interest of efficiency it 
is sometimes desirable to program that trajectory which requires the least 
amount of energy. This paper considers the minimum energy movement for 
a spring-like actuator abstracted from muscle mechanics and energetics. It 
is proved that for this actuator a bang-coast-bang actuation pattern 
minimizes the energy expenditure. For some parameter values this pattern 
is modified by a singular arc at the first switching point. A surprising 
limitation on the duration of coast is demonstrated. Some relaxations of the 
restrictions underlying the spring model are shown to preserve the bang- 
coast-bang solution. 
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Though considerable effort has been expended in the study of the 
human motor system, the execution of even simple movements is not well 
understood. One current theory holds that movements are memorized in terms 
of final position CBizzi et al.I. The organism selects length-tension 
curves of agonist and antagonist by adjusting the innervation levels so 

that the intersection of those curves occurs at the desired position 
[Fe I dman]. 

This process is illustrated by the hypothetical length-tension 
curves of agonist and antagonist muscles in figure 1. Suppose the system 
is currently at length Lq under innervation rates g| for the agonist and nj 
for the antagonist. If the innervation rate of the agonist is changed to 
g 3 , a different agonist length-tension curve is selected and the 
equilibrium length shifts to Lj. Assuming no delay in tension development 
and ignoring velocity effects, the arrow in the figure indicates the 
tension course. There is an isometric buildup of tension from P 0 to P 2 
followed by an isotonic decay to P|, where the tension in agonist balances 
the tension in antagonist. The final position theory maintains that the 
position L| can be reached independent of starting position merely by 
setting rates nj and gs. This theory is Interesting from a manipulation 
viewpoint because it obviates the need for precise trajectory calculation. 

There are many choices of agonist-antagonist length-tension curve 
pairs that have Lj as equilibrium position. One choice that could be 
expected to require less energy is n 2 and g 2 , which minimize the isometric 
tensions. More generally, it is conceivable that some complex sequence of 
innervation rates Inj,g } l might require less energy than a scheme which 






Figure_JL. The equilibrium point of the intersecting 

length-tension curves of agonist (g labels) and 

antagonist (n labels)' shifts from L Q to when the 

firing rate of the agonist is raised from q to q 

^1 ^ 3 

and the antagonist rate remains at n^. 


Tension 



Figure 2 . Length-tension curves from the cat soleus 
muscle (Rack and Westbury). 




selects the final length-tension curves immediately. The determinaton of 
this optimal innervation pattern is the focus of this paper. 

I. Muscle Mechanics 

The investigated properties of muscle present a too complicated 
view for analytic treatment. The plan in my research has been to simplify 
the muscle mechanics until an analytic solution to the optimal energy 
problem is possible, then to examine if the nature of the solution is 
changed by adding some of the excluded muscle properties. A full-blown 
formulation of the problem Mould require numerical methods for solution, 
and I intend eventually to carry out this analysis. 

The length-tension curves in figure 1 are themselves a 

simplification of actual length-tension curves such as in figure 2. The 

* 

simplification arises from extracting only the linear portions at short 
lengths. There is some question as to what portion of the length-tension 
curves are used in actual movement. Some authors [Zierler, Hill 1970, Cook 
and Stark] maintain that the active portions occur near L 0 , the length at 
uhich there is maximum isometric tension and uhich Zierler also calls the 
rest or natural length of muscle. Collins et al. 11975] on the other hand 
report that the linear portions are used in eye movement. 

For the present we accept the hypothetical length-tension curves of 
figure 1. The curves are also assumed parallel, as reported by Collins et 
al.; Rack and Uestbury, houever, report a decrease in slope with firing 
rate for this linear portion (figure 2). The curves in figure 1 lead to a 
model of muscle as a spring with variable zero setting. The slope K of the 









PAGE 5 


curves represents the spring constant, and the variable zero setting L t 
corresponds to the selection of firing rate. The force exerted by a muscle 
is thus K(L-L,). 

An important simplification of muscle properties is to exclude the 
series elastic component. Ue also exclude the parallel elastic component 
and the active damping. The resulting muscle model is presented in figure 
3. The equation of motion for the spring system of figure 3 let 

m d 2 x ■ - b dx + kg(x.-x) - k n (x-x n ) (1) 

dt 2 « 

t 

Define a control variable X and a state variable U as belou. 

. t x, ] [ x ] [ u, ] [ x. 1 

X - [ ]-[} U-t ]-[ ) 

t x 2 ] til t u 2 ] C x„ ] 

Setting the mass m » 1, the state variable representation of the spring 
system is: 

. t 0 1 1 [00] 

X - t ] X + I ] U (2) 

t -k g -k n -b ] t k g k„ 1 

More compactly, 

X - A X + B U (3) 

- f(X,U,t) 

II. Muscle Energetics 

The energy E expended during movement equals work plus heat. The 
work U may be subdivided into conservative uork performed on the mass m and 





Simplified muscle model with equation of motion: 

• • 

m x = - b x + k (x -x) - k (x-x ) 

9 g n n 

b * coefficient of passive damping 
b = spring constant of the agonist 

spring constant of the antagonist 
m = mass 
x * position 

X g = var ^ a ^ le zero setting for the agonist 
x n “ var ^ a ^^ e zero setting for the antagonist 
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nonconservative work performed on the viscous element. The isometric heat 
Qj is given off in maintaining the muscle at a particular tension Pg. The 
rate of energy expenditure is thus 

E - P 0 v (power) + a P 0 (maintenance heat rate) 

■ (x 2 + a) (k g (uj-xj) + k n (x|-u 2 )) (4) 

where v is velocity and a is the maintenance heat coefficient# The two 
force terms have been summed because each contributes to energy toes# We 
have excluded the shortening heat because the active damping was also 
excluded, and because there may be a theoretical relationship betueen the 
two [Huxley, Caplan), The transient characteristics of heat production 
have also been excluded. 



III. The Euler-Lagrange Equations 

The task now is to find the time varying control U(t) that 

minimizes the energy used in moving between two points in a fixed interval 

of time. Let V represent the energy consumed in applying the control U to 

yield the trajectory X. The problem of minimizing V is readily approached 

by techniques of modern control theory. The fundamental equations that the 

optimal control U(t) must satisfy are derived from a theorem from the 

calculus of variations . This theorem states that in o^der to find the n- 

vector X(t) that minimizes V(X), where 

r*f 

V(X) - / L(X,X,t) dt (5) 

J*. 

subject to the constraint relations 

g- ( (X, X, t) - 0 I » 1,... ,m < n (6) 
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then X(t) satisfies the Euler equations 

3L* (X.X.t) - d 3L’(X.X.t) « 8 1-1.n (7) 

3xj dt 9>tj 

where 

LMX.X.t) - L(X,X,t) + 2 Xj(t)g,(X,X,t) (8) 

and Xj(t)'e are the multiplier functions [Schultz and lie I sal. 

Applying this theorem to the optimal control problem, the state 
* 

equations X - f(X,U, t) represent the equality constraints. L represents 
the rate of change of energy E. The Hami I ton!an H - L +• X T f repreeents 
(8), where X^ - IXj X2) • By applying the Euler equation first for X and 
then for U, it can be shoun that the minimizing U(t) satisfies the 
following two Euler-Lagrange equations [Schultz and Me I sal. 

X - -H x (9) 

H u - 8 (10) 

IV. The Minimum Principle 

Because L is linear in the control U, there will not generally 
exist a minimum energy solution. To obtain a realistic solution, 
constraints must be placed on the control. The solution in this case will 
lie on the constraint boundaries [Bryson and Hoi. Constraints on U, 
however, make it impossible to differentiate H with respect to U. 

The minimum principle of Pontryagin makes it possible to proceed 
from this point. Pontryagin showed that even if the control is 
constrained, one still obtains a minimal solution by finding the u° « 
u°(X,X,t) to minimize the Hamiltonian H, but by inspection rather than by 
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differentiation. After finding the minimizing u°, one forms H° • 
H(X,u°,X,t) and then solves the following two equations [Schultz and 
Me I sal. 


3H° 

(11) 

ax 


- 3H° 

(12) 


ax 


There are two natural constraints that fall on the control U. 

First, the spring cannot push. 

U, - x, £ 0 (13) 

*1 - u 2 > 0 (14) 

Second, springs have a maximum tension that they can exert. Uithout this 
constraint the solution would involve an infinite impulse. For the moment 
we assume the maximum tension is constant and independent of length: 

Uj - X| s. C| (15) 

xj - u 2 £ c 2 (16) 


where C| and c 2 are constants. The case of maximum tension varying with 
length is deferred until section XI. 


V. A Bang-Coast-Bang Solution 

To facilitate inspection of the Hamiltonian, we expand H - X T f + L 
into three lines, the first depending on u t , the second on u 2 , and the 
third on neither control. 

H ■ kg(uj-xj) (tt+x 2 +X 2 ) (17) 

+ k n (xj-u 2 ) (o+x 2 -X 2 ) 

+ x 2 (Xj-bX 2 ) 
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To minimize H with respect to Uj, ue observe that if a+X 2 +X 2 > 0 then H is 
minimized when U|«X|. If a+x 2 +X 2 < 0 then H is minimized uith U|-X|+c t . 
Similarly, it can be shown for u 2 that when a+x 2 -X 2 < 0 the minimizing u 2 
lies at X|-c 2 ; otherwise u 2 is at xj. Combining these results, one finds a 
bang-coast-bang solution to the minimum energy for muscle movement. 

Case 1: X 2 < -Io+x 2 ) 

Then Uj-Xj+Cj, u 2 -Xj. 

Case 2: IX 2 I < (a+x 2 ) 

Then Uj-xj, u 2 «xj. 

Case 3: X 2 > («+x 2 ) 

Then u ( -X|, u 2 «x t -c 2 . 

% 

VI. The Solution Equations 

f 

Substituting the minimizing u° into H, one obtains three functions 
corresponding to the three cases. 

Case li ■ k|C| (ot+x 2 +X 2 ) + x 2 (Xj-bX 2 ) 

Case 2s H° - x 2 (X|-bX 2 ) 

Case 3: H° - k„c 2 («+x 2 -X 2 ) + x 2 (X|-bX 2 > 

The differential equation (11) and its solution becomes for the three 
cases: 
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,/-b(t-t 0 ) . ,, -b(t-t 0 ), ,, B4 

x 2 (t) - x 2 (t 0 )e u + kcCl - e w ) (18) 

b 

1 ♦ 

x,(t) - X|(t 0 ) + 1 (x 2 (t(j) - x 2 (t)) + kc (t-t 0 ) (19) 

b b 

where c-cj and k-k g for case 1* c-0 for case 2; and c—c 2 and k-k„ for 
case 3. The time t 0 represents the starting time. The differential 
equation (12) and its solution ares 

(80) t 0 ] 

A - I 1 X - ( 1 

t -1 b ) ( kc 1 

A|(t) - X,(t 0 ) (20) 

X,(t) - X 2 (t 0 )e b(t " to) 4- Ai (tn)+kc (1 - e b(t “ to) ) (21) 

b 

_ where c and k have the same meaning as above except c-+c 2 for case 3. 

Since XiCt) is constant, it appears henceforth as A] without a time 
dependence. 

VII. The Extremal Versus Singular Solution 

It is proven in appendix A that there are exactly three events in 
the extremal bang-coast-bang solutions an acceleration period, a coast 
period, and a deceleration period. No other combination of bangs and 
coasts is minimizing. However, a nonextremal minimizing solution may arise 
from a singular arc at the switching points. The Hamiltonian (17) has the 
curious property that if X 2 "|a+x 2 | then the corresponding control may take 
on any value and still minimize H. If a control can be found to maintain 
X 2 *|a+x 2 | for a finite time interval, then a non-extremal solution to 






energy minimization might exist. This situation is called a singular arc 
and arises from a performance index linear in control but quadratic in 
state [Bryson and Hoi. To maintain X 2 B la+x 2 | for a finite time interval, 
all time derivatives of the tuo suitching curves must be zero: 

d"tX,+g+x,) - 0 n*0 (22) 

dt" 

d"(X,-a-x,) - 0 n*0 (23) 

dt n 

Carrying through the analysis for a singular arc at the first suitching 
point (22), the time varying force during the singular arc is: 

k,(U|-X|) - Xj + b« + 3bx 2 (t,)e 2b(t ' t,) + 3(ba+X,) (e 2blt “ t|> -l) (24) 

Unfortunately no sufficient condition has yet been developed to 
test uhether a singular arc is minimizing, and one must compare values of 
the performance index for specific parameter values for the singular arc 
solution versus the extremal solution. Depending on the choice of X| the 
force (24) takes one of the three forms in figure 4. 

Of these forms only 4C has been found minimizing for some parameter 
combinations. To search for such combinations, a set of parameters uas 
initially deduced from Rack and Uestbury (table 1). The elapsed distance 
and the elapsed time tf are variable and have been chosen as 0.2 cm and 
0.4 sec respectively. The initial and final velocities are assumed zero. 
For the extremal solution (bang,coast,bang) there result 8 nonlinear 
equations in 8 unknowns from (18)-(21) and the initial conditions lx|,tf). 
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Table I 

k ■ 2 kg/cm 

in - 0.2 kg (plausible value) 

b/m * 3.16/sec (chosen to give {’■0.5) 

c ■ 1 cm 

a - 0.1 cm/sec (deduced from Mo I edge) 

* 

For the singular solution (bang,4C,coast,bang) 15 nonlinear equations in 15 
unknowns result from (18) — (24) and the initial conditions. The equations 
were solved numerically by Neuton-Raphson and gradient methods. Individual 
parameters were varied and energies of movement computed from (5). Solving 
(5), the energy for the extremal solution is 


E ■ kcOxj (t|)+dx|(t 3 )) + kcoOt 1 +3t 3 ) (25) 


where t| is the switching time from acceleration to coast, t 3 is the time 
at the end of deceleration, dt| is the duration of acceleration, dx|(t|) is 
the distance moved during acceleration, and dx t (t 3 ) is the distance moved 
during deceleration. For the singular solution, the energy is 


(26) 

where tj is the switching time from acceleration to the singular arc 4C, t 2 
is the switching time from 4C to coast, and t 4 is the time at the end of 
deceleration. The force k g (uj-xj) is given by (24), while the velocity 
x 2 (t) is 


E ■ kc(dxj(t|)+dx|(t*)) + kca(dt|+dt 


r 

4) + / k,(U|-X|) (x 


2 (t)+a) dt 


x 2 (t) - x 2 (t|)e 


2b(t-t,) 


2b(t-ti) 

ba+\, (e 1 -1) 

2b 


(27) 


The energies for the extremal versus the singular solution are 
compared in tables 11 a-g; the units are kg cm/kg wt. In table lie 
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Table Ha 


k 

Singular 

Extremal 

7 

impossible 

0.708 

8 

0.66117 

0.66106 

*10 

0.6173 

0.6164 

20 

0.5640 

0.5627 

30 

0.5513 

0.5509 

40 

0.54554 

0.54552 

50 

0.5423 

0.5426 

100 

0.536 

0.537 


Table IIc(b*fa 


k 

Singular 

Extremal 

7 

impossible 

0.664 

8 

0.632026 

0.632023 

*10 

0.6173 

0.6164 

14 

0.6442 

0.6440 

20 

0.708 

0.725 


Table 

lie 



Singular 

Extremal 

0.35 

impossible 

0.814 

0.36 

0.76123 

0.76122 

*0.4 

0.6173 

0.6164 

0.45 

0.512 

0.511 

0.5 

0.444 

0.445 

0.6 

0.358 

0.369 


Table lib 


b 

Singular 

impossible 

Extremal 

2.3 

0.545 

2.4 

0.552036 

0.552034 

2.6 

0.5675 

0.5674 

3.0 

0.6020 

0.6014 

*10 

0.6173 

0.6164 

4.0 

0.7045 

0.7030 

5.0 

0.820 

0.824 

10.0 

1.45 

1.51 

15.0 

2.13248 

2.13246 

15.557 

impossib1e 

2.4 


Table lid 


C 

Singular 

Extremal 

WE * vj % 

impossible 

0.708 

0.8 

0.6612 

0.6611 

*1.0 

0.6173 

0.6164 

2.0 

0.5640 

0.5627 

3.0 

0.5513 

0.5508 

4.0 

0.54554 

0.54551 

5.0 

0.5427 

0.5437 

6.0 

0.5401 

0.5407 

Table 

IIf(ksl6. b=4) 

• 


Xf 

Singular 

Extremal 

0.006 

impossible 

0.0041 

0.05 

0.0632 

0.0634 

0.1 

0.1925 

0.1942 

*0.2 

0.6647 

0.6681 

0.3 

1.4636 

1.4620 

0.35 

2.0120 

2.0072 

0.4 

2.6878 

2.6825 

0.54 

impossible 

6.4598 
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Table IIg (k=16. b*4) 


a 

Singular 

Extremal 

0.06 

0.6112 

0.6161 

*0.1 

0.6647 

0.6681 

0.2 

0.7974 

0.7981 

0.25 

0.8633 

0.8631 

0.5 

1.1903 

1.1882 

1.0 

1.8393 

1.8382 

1.7 

2.748417 

2.748411 

1.8 

impossible 

2.8784 


the parameters k and b are varied simultaneously but at a fixed damping 
ratio of 0.5. In tables Ilf-g the parameters k and b are respectively set 
at 16 and 4 rather than at the table I values where the extremal solution 
is minimizing over the whole range of x f and a. The initial values from 
table I are starred in table II. A singular solution becomes minimizing 
with high values of k, b, c, and tf, and with low values of a and x f . As 
the parameters cause the coast time to approach zero (higher b and x f , 
lower k, c, and t f ), the singular and extremal solutions become identical 
because the 4C portion vanishes. 

For the extremal solution it is proved in appendix B that there is 
an upper limit on the duration of coast. It is tempting to speculate that 
for longer coast durations a singular solution becomes minimizing, but the 
singular solution in table II is not always minimizing under these 
conditions. Perhaps a different combination of bangs, coasts, and singular 
arcs would then be minimizing, but this remains an open question. Some 
combinations can be proved impossible, such as (bang, 4C, coast, 4A-C, 
bang). 





VIII. Spring Model Relaxations: X Dependencies 

A natural question is whether the minimum energy solution is 
changed by incorporating a more realistic muscle model. For those 
relaxations of the spring model involving only X dependencies, the answer 
18 that the solution remains bang-coast-bang. The reason is that the 
Hamiltonian H remains linear in the control U, and the minimization of H 
with respect to U occurs at fixed X. Mhether the solution also remains 
acceleration-coast-deceleration needs to be determined for each case. 

Relaxations of the spring model involving X dependencies include 
the following. 

1. Position Limits on Tension 

For real muscle the maximum isometric tension varies with position 
(figure 2). This makes Cj and eg into functions of xj, but the controls 
will still fall at the extremes wherever they are. 

2. Velocity Limits on Tension 

Actual muscle exhibits a hyperbolic force-velocity relation. If p 0 
is the isometric tension, then the maximum force P that can be produced for 
a velocity v is [Hill 19381 (see figure 5): 

P - P 0 - v (Pn+a) (28) 

v+b’ 

The term (Po+a)/(v+b*) can be considered the coefficient of active damping. 
The coefficient a has been determined as .25 P 0 j the force P then becomes 
Po“l*25P 0 v/(v+b’). The literature conflicts on the value of active damping 
during lengthening. For consistency with the shortening heat (below), it 




Tension 


.?A9 ur e 5 . Tension dependence on velocity (Hill 1938) 



length 


ll^9 ure 6* Hypothetical length-tension curves with 
the property that at any given length the slopes are 
the same for all choices of u . 
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is assumed the same as active damping during shortening. 

Associated with the active damping is an extra heat expenditure 
above the isometric heat due to shortening. This shortening heat rate is 
[Hill 1964]: 

Q 8 - .16 P 0 v + .18 Pv (29) 

The isometric heat rate Qj remains «Po, but the power is now Pv. 
Substituting the sum of 9pring forces for the isometric tension P 0 and (28) 

for P, the energy rate i sc 

« • 

L - (k,(U|-X|) + k n (u 2 -x,))(« - 8.09><2 + 1.45b» ) (30) 

X2+b* 

Similarly it can be shou that the equation of motion ist 

x 2 » - b x 2 + (1 - 1.25xj>) (k f (ui-xi) - k n (x|-u 2 )) (31) 

x 2 +b’ 

When these terms are combined to form the Hamiltonian, the control is seen 
to remain linear. Hence the solution is once again bang-coast-bang. 

3. Soring Constant Variations with Position 

One way of bringing the simplified length-tension curves of figure 
1 closer to those of figure 2 is i I lustrated in figure 6. The spring 
constant k n varies with position, but at any given position the constant k„ 
is the same for all controls u 2 . Under these conditions the solution 
remains bang-coast-bang. 

4. Parallel and Series Elastic Elements 

The incorporation of these elements into the model is depicted in 
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figure 7. Since the parallel elastic element depends only on position, it 
does not change the solution. The series elastic elements and the active 
springs may be replaced uith equivalent springs uith constants 
K f *«k g k 8 /(kg+k,) and k n ’-k n k,/(k B +k # ). This modification also has no effect 
on the solution. 

IX. Spring Model Relaxations: V Dependencies 

In figure 2 the spring constant k n is seen to vary uith firing rate 
at any fixed position. The linear portions of these length-tension curves 
uhen extended seem to intersect at a common point (figure 8). In this 
circumstance the spring constant k„ is Tq/(u 2 -x 0 ). The u z terms of the 
Hamiltonian H become: 

H* * To (xi~Uj) (fl+xj'Xjl 

Uj-Xq 

■ To(a+x z -X2) ( xi-Xft— 1) (32) 

U2“Xq 

If a+* z -X2 > ® then H’ is minimized at U 2 "*i: otheruise U 2 »X|-c z . That is 
to say, the solution for u z is exactly the same as in section V. A similar 
analysis holds for uj. Thus the minimizing pattern is also bang-coast- 
bang. 


F igure —7_. An expanded muscle model incorporating series 
elastic elements k and a parallel elastic element k . 



£igurje_8. The length-tension curves of Rack and Westbury 
when extended meet at a point. The dependence of k 

n 

can be characterized by T /(u -x ). 


on u 



Appendix A 

In this appendix it Is shown there are only two bangs and one coast 
in the extremal solution: one acceleration, followed by one coast period, 
terminated by one deceleration. No other combinations of coasts and bangs 
are possible. To demonstrate this is the only possible combination, it Is 
necessary to examine the switching curves and their time derivatives. 

The first lemma shows that once the control has passed from 
acceleration to coast, then the control cannot return to another 
acceleration but must proceed to deceleration. The second lemma shous that 
once deceleration has started, the deceleration must continue until the end 
of the movement. This proves that the acceleration-coast-deceleration 
combination is the only possible one. In the following it is presumed that 
the movement starts with acceleration in the positive xj direction. Hence 
all velocities are positive. 

Lemma It After acceleration, the glide period cannot double back to 
another acceleration. 

Proof: The proof of this lemma proceeds by examining the time derivative 
of the acceleration-coast switching curve (henceforth referred to as the 
elope of the switching curve). The slope of this curve is initially 
positive at the transition from acceleration to coast. In order for 
another acceleration to follou the coast period, this slope must become 
negative, leading to a contradiction. 

At the first switching time t| the acceleration-coast switching 
curve Is zero. 
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^ 2 ^ 1 ^ + + *0 

After the acceleration period, the coast equations are: 


(Al) 


* in * i b(t-tj) . btt-tj). 

^(t) ■ X 2 11 1 )e + 3^J_C 1 - e ) 

b 


x 2 (t) ■ - (Aa(tj)+«)e 


-blt-t|) 


(A2) 


(A3) 


Thue 


X2^)+®+><2^) * X2<t|)e * + Xid-e * ) 

b 

..... . -b(t-tj) 

+ # ■ (X2(tj)+«)e 


(A4) 


The elope of this switching curve is: 

e b(t ' t,, (bX 2 (t,)-X,) + b(X 2 (t,)+a)e" b(t " t,) (A5) 

# 

At t«tj, the 81 ope of the switching curve Is: 

2bX 2 (t|)-Xi+ba > 8 (AG) 

One can show this quantity cannot be less than zero. Next, suppose the 
coast doubles back to another acceleration. At some point the slope must 
go through zero. This time t is found from (A5) as: 

2b 

e 1 - b(X»(t,)+a) > 1 (A7) 

X)-bX2(11) 

Case 1: X|-bX 2 (t|) > 0. 

Then X2(t|)+a>0, contradicting (Al). 

Case_2j_ X|—bX 2 ^tj) <0. 

Crossmultiplying (A7) and collecting terms, 

2bX2(t|)-X|+ba<0, contradicting (A6). 
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Thus after acceleration, the coast period must eventually arrive at the 
deceleration suitching point. 

Lemma 2: The movement is locked in deceleration until the end. 

Prooft It Mill be shown that if decelsration ever switches to coast, then 
the slope of the coast-deceleration suitching curve requires an immediate 
return to deceleration. Hence the movement is locked in deceleration until 
the end. 

Suppose there is a time t 3 when deceleration suitches to coast. At 
this point the coast-deceleration suitching curve is zero. 

Jijitji-e-xjltj) ■ 0 (A8) 


The coast suitching curve X 2 (t)-«-x 2 (t) is 


\ , b(t-t 3 ) ,. b(t-t 3 ) . 

X 2 ^ 3 )e + ) - a - x 2 (t 3 )e 

b 


-b(t-t 3 ) 


(A9) 


The 81 ope of (AS) ist 


b(t-t 3 ), ...... . . -b(t-t 3 ) 

e 3 (bX 2 (t 3 )—Xj) + bx 2 (t 3 )e 3 


(A10) 


At time t 3 the slope (A10) is bX 2 (t 3 ) —Xt**’bx 2 (t 3 ). This is positive since 
X2(t 3 )>0, x 2 (t 3 )>0, and X )<0 (lemma 3). This means that deceleration uould 
bounce off the coast boundary and immediately continue the deceleration. 
Furthermore, since the slope is positive, the deceleration uould not 
immediately suitch back to coasting, causing chattering. 






Lemma 3: X] < 0. 

Proof: At the second snitching point t 2 we have 


From (A2) and (A3), this becomes 


A 2 (t,)e b(t2 ‘ t,, + i LL (l-e b(t2 “ t,) ) - a + (A 2 (t,)+a)e’ b(t2 " t,) 

b 


Rearranging, 

\ ii h(t 2 -t|) . . . . b(t 2 -tj) -b(t 2 -t|). M -b(t 2 -t 

XjJl-e ) - - X 2 (tj)(e * 1 +e * 1 ) + «(l-e 

b 

From (Al) and (A3) ue find an expression for X 2 (t). 

X 2 (t,) - - « - k i c 1 (l-e‘ b(trt ° ) ) 

b 

Substituting into (A13), 


i l (l-e b(,2 " t|, | - k-cU-.'” 11 ''* 01 ) (« b(,i ' t|1 t«~ l ’ l<2 ' t|1 )) 


Thus 




1-8 


b(t 2 -tj) 


Since the numerator is positive and the denominator is negative, X 
negative. 


Taken together, these lemmas show that acceleration.passes through 
deceleration. There is no possible variation in this scheme. It 
possible to shou the movement cannot start by coasting folloued by 
acceleration. 


PAGE 25 


(All) 

0 (A12) 

,} ) (A13) 

(A14) 

* 

(A15) 


(A16) 

is 

coast to 
s also 
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Appendix B 

A surprising limitation on the value of t 2 “tj, the duration of the 
coasting time, has been found. The suitching curve during acceleration is: 


X2(t)+tt+x 2 (t) - X2 1 1 0 ) e* 3 * * Ai+k t ci (l-e b ^ * 0 *) + a 

b 


-b(t-t 0 ). 
+ k.ci(l-e w ) 


The slope of this switching curve is: 


(Bl) 


b(t-t 0 ) % . -bCt-t 0 ) 

e u (bXjttol-Xr^fCi 1 + kfCje 


CB2) 


At the first suitching time t|» the suitching function (Bl) is zero. 


Rearranging (Bl) for t«t|, 

e b(t|-t°) (X9 ( tn) . ( x, +k e ,) ) . . ( X.-^kc) - 

b b 

Substituting (B3) into (B2), the slope at t| is: 


- k-Cj_(l-e~ b ^ , ” t °S (B3) 

b 


- 2k,c,(l-e" b(trto) ) - Xi - ba 


(B4) 


Substituting for X| from (A16), 


bltj-tj) 


(B5) 


Since slope(ti)>0 and since the denominator is positive, so is the 
numerator. 


M 1 aV b, * r * 0 W' b,,r * l ’-. bt, *-‘> , ) + 2b„ i e 


(B6) 
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_ _ -b(t 2 -t.) b(t 2 -t.) , 

Because 2+e -e Is a decreasing function of t 2 , at some 

point CBS) becomes zero. Solving then for e b * 2 


bo 


+ V(bp+k t ci(l-e *°*)) 2 + (k t c (l-e to *)) 2 ' + 1 


(B7) 




As t 2 increases, tj Mill decrease. However, t| does not decrease enough to 
offset the effect of the t 2 increase. If p-0, (B7) reduces to 


Strangely, in this circumstance t 2 -tj depends only on b. 


(B7a) 
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